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We study the distribution of residence time or equivalently that of "mean magnetization" for a family of Gaussian Markov 
processes indexed by a positive parameter a. The persistence exponent for these processes is simply given hy 9 = a but 
the residence time distribution is nontrivial. The shape of this distribution undergoes a qualitative change as 9 increases, 
indicating a sharp change in the ergodic properties of the process. We develop two alternate methods to calculate exactly but 
recursively the moments of the distribution for arbitrary a. For some special values of a, we obtain closed form expressions of 
the distribution function. 

that includes the diffusion equation [pi. 

Recently it was argued |^,^ that given this stochas- 
tic process y{t), it might be useful to investigate a more 
general quantity, namely the "residence time distribu- 
tion" , whose limiting behaviour determines the persis- 
tence exponent. This is the distribution f{r,t) of the 
random variable, r{t) = 7/0 9[y{t')\dt' where 9[x) is the 
Heaviside theta function. Thus r(t) is simply the frac- 
tion of time spent by the process y{t) within time t on 
one side of zero. It was shown in ref. [E5| that for any 
Gaussian stochastic process, the distribution f{r) is in- 
dependent of time t. For Gaussian processes with zero 
mean, the symmetry r <-> (1 — r) indicates that the func- 
tion /(r) is symmetric around r = 1/2. Also in the limit 
r -^ (and symmetrically for r ^ 1), the function /(r) 
is clearly the probability that the process remains only 
on one side of zero and hence is proportional to per- 
sistence. This indicates that as ?' -^ 0, the function 
/(r) must behave as ^ r^^^ (and as ^ (1 — r)^~^ for 
r -^ 1), so that f{r)dr ~ t~^ as r -^ or 1. A somewhat 
more convenient variable is the "mean magnetization" 
pjl, m,{t) — 2r{t) — 1, whose range is [—1, 1] and whose 
distribution function, P{m) = ^/[(l + "^)/2] is symmet- 
ric around to = and behaves as P{m) ~ (1 ± m)^^^ 
near tti = ±1. 

The distribution P{m) is known exactly for the process 
that represents the position of a one dimensional Brown- 
ian walker [ p6[ . Lamperti p7| derived an exact expression 
of P{m) for a class of renewal processes where successive 
zero crossing intervals are statistically independent. Re- 
cently a special case of Lamperti's results [pTJ, when the 
successive intervals are distributed according to a Levy 
law, was rederived by Baldassari et. al. ||2^ by a differ- 
ent method. The distribution P{m) has been determined 
numerically for diffusion equation |25| ] and for interface 
growth models [£1|. Besides, moments of P(m) have been 
determined analytically for diffusion equation under the 
independent interval approximation P4| . 

The distribution function P{m) provides a somewhat 
more detailed information on the statistical nature of 



I. INTRODUCTION 

The problem of persistence in spatially extended 
nonequilibrium systems has recently generated a lot of in- 
terest both theoretically ||l|-Rl and experimentally [p|-p^ . 
These systems include the Ising or Potts model undergo- 
ing phase ordering dynamics |lj-|j,|ll|-|l6| , simple diffusion 
equation with random initial conditions ||5|,|6|, several re- 
action diffusion systems in both pure (l^ and disordered 
|l8| environments, fluctuating interfaces [[l9|-pi[, Lotka- 
Volterra models fE2| and inelastic collapse of a randomly 
forced particle ||23| . In many of these systems the spatial 
degrees of freedom of the original many body problem can 
be integrated out and the problem of persistence effec- 
tively reduces to the calculation of the probability Pq (t) 
of no zero crossing upto some time t of an effective single 
site stochastic process y{t). 

In most cases of interests, this probability decays as a 
power law for large time, Po{t) ~ t~ , where the persis- 
tence exponent is nontrivial. This nontriviality can be 
traced back to the fact that once the spatial degrees of 
freedom are integrated out, the effective single site pro- 
cess y{t) becomes non-Markovian. For a non-Mar kovian 
process, it is well known that the calculation of any his- 
tory dependent quantity such as persistence (no zero 
crossing probability) is extremely hard [ p9pC| ]. As an 
example, for the diffusion equation with random initial 
condition, the effective single site process y{t) is a Gaus- 
sian non-Markovian process characterized by its two-time 
correlator, {y{ti)y{t2)) = [4*1*2/(^1 +12^]'^^'^ where d is 
the spatial dimension Q. Even for this simple case, the 
corresponding persistence exponent 6 is nontrivial and is 
known only numerically and approximately by analytical 
methods ||,g^ but not exactly so far. Though recently, 
an exact series expansion result for the exponent 9 has 
been derived for arbitrary smooth Gaussian processes 



the stochastic process y{t). For example, in the con- 
text of diffusion equation it was pointed out by New- 
man and Toroczkai pq ] that an interesting information 
can be extracted from the shape of the function P{rn). 
For diffusion equation, the exponent 9(d) (which con- 
trols the shape of the function P{m) near m — ±1) in- 
creases monotonically with space dimension d. There 
exists a critical dimension dc where 9{dc) = 1 such that 
for d < dc, 9 < I and the function P{m) diverges as 
771 —> ±1, has a minimum at ttt, = and is concave up- 
wards in the range [—1, 1]. On the other hand, for d > d^, 
6* > 1, the function P{m) goes to zero as 777 — > ±1, has a 
maximum at 777 = and is convex upwards in [— 1 , 1] . The 
peak of the distribution shifts from the edges 771 = ±1 to 
the center m = a.s d increases through dc- Thus for 
d < dc, the most probable configurations of the process 
y(t) are the ones which do not cross zero whereas such 
configurations are least probable for d > dc, signalling 
the existence of a sharp change in the ergodic proper- 
ties of the diffusion field. Such a detailed information is 
not contained in the persistence exponent 9. In rcf. P5[ , 
dc for diffusion equation was approximately determined, 
dc ~ 36. 

These useful informations contained in P{m) of the 
diffusion equation have so far not been possible to derive 
exactly mainly due to the non-Markovian nature of the 
single site Gaussian process. It would therefore be useful 
to find and study some simpler Markovian Gaussian pro- 
cesses with some tunable parameter (which would play 
the similar role as the spatial dimension d does in dif- 
fusion equation) where exact calculations can be per- 
formed. In this paper we study the magnetization distri- 
bution P{m) of a family of such Gaussian Markov pro- 
cesses parametrized by an index a. By varying this pa- 
rameter a, the persistence exponent 9 for this process 
can be varied continuously. The Markovian nature of 
the process also makes many exact calculations possible. 

The Markov process y{t) that we study in this paper, 
satisfies the following stochastic Langevin equation. 



dt 



(1) 



where 'q{t) is a Gaussian white noise with {rj{t)) = and 
{r]{ti)ri{t2)) = 5{ti — ^2) and a is a positive parameter. 
There are various physical processes that are described by 
the above Langevin equation. For example, for a — 1/2, 
y{t) represents the position of a one dimensional Brown- 
ian random walker. For a — 1/4, y{t) can be interpreted 
[ pT[ to be the "total magnetization" of a Glauber chain 
undergoing zero temperature coarsening dynamics after 
being quenched rapidly from infinite temperature. 

The persistence of the process y{t) is simply the prob- 
ability for this process not to cross zero upto time t 
and decays as t^^ for large t. The exponent 9 for this 
process can be trivially computed, 9 = a. The sim- 
plest way to derive this is to define a new time vari- 



able t' — t^" such that the equation of motion becomes 
dy/dt' — ({t') where the new noise C{t') has zero mean 
and (C(ii)C(4)> = ^it'i - ^2)- But this is simply the 
equation of motion of a one dimensional Brownian walker 
whose probability of no return to zero upto time t' de- 
cays as ^ l/Vt^ ^ t^"- Thus the persistence exponent 
for y(t) is simply 9 — a. 

However we show in the rest of the paper that even 
though the persistence exponent 9 = a trivially for this 
process, the magnetization distribution Pirn) is non- 
trivial. In fact, as a and hence 9 is increased, the 
shape of P{m) changes from concave upwards to con- 
vex upwards. For a — 1/2 (i.e., for ordinary Brownian 
walker), the distribution P{m) was already known ex- 
actly, P[m) = l/(7r-\/l — 777^) ll^. For general a, while 
we have not been able to determine the full distribution 
P{m) in closed form, we demonstrate below by two com- 
pletely different methods that the moments of P{m) can 
be calculated exactly. In the first method we generalize 
the formalism developed by Kac |3^ for a = 1/2 case to 
arbitrary a. In the second method, we use the formalism 
recently developed in the context of diffusion equation by 
Dornic and Godrechc 124] using independent interval ap- 
proximation (IIA) . Wc point out however that while this 
latter method yields only approximate results for diffu- 
sion equation |Q, it gives exact results for the Markov 
processes that we study in this paper. 

The paper is organized as follows. In section (II), we 
generalize Kac's formalism for a — 1/2 to arbitrary a 
and derive an exact recursion relations satisfied by the 
moments of P{m). In section (III), we rederive the same 
results by using an alternate IIA formalism. In section 
(IV) we use the formalism developed in sections-II and 
III to obtain explicit results for the distribution of mean 
magnetization for some special values of the parameter 
a. Finally we conclude with a summary and discuss the 
relative merits of the two formalisms and some applica- 
tions. 



II. METHOD I: GENERALIZATION OF KAC'S 
FORMALISM 

We consider the Gaussian process y{t) evolving 
stochastically via Eq. |l] and define the "mean magnetiza- 
tion", m{t) — J Jq V[y{t')]dt' , where the functional V{y) 
in our case is simply, V{y) = sgn{y). Let G{y,t \ y',t') 
denote the propagator of the process, i.e., the probability 
that the process takes the value y at time t given that it 
was at y' at time t' < t. This can be easily computed for 
our process and is given by. 



Giy,t\y',t')^ 



^{y-y'fmt^ 



^2lT{f^^ 



t'2") 



(2) 



Following Kac |31|, we define the moment generating 



function, 



Etf- 



n=0 



where f„ are the moments defined by: 



^n = ((/ V[yit')]dtr)- 
Jo 



(3) 



(4) 



To compute the moments Vn, it is useful to first define a 
set of functions Q„(j/,i) via the recursion relation, 

Qn+i{y,t) ^ dt' dy'G{y,t \ y' ,t')V{y')Q^{y\t') 

Jo J -00 

go(2/,i)=G(y,t|0,0). (5) 

It can then easily be checked that, 

/oo 
Qniy,t)dy. (6) 

-00 

Using Eq. || and Eq. ^, we finally get, 

dyQu{y,t) (7) 

where Qu{y^t) is the generating function, 

00 



Thus the moments of the mean magnetization m can 
be computed exactly from Eq. M provided we can evalu- 
ate the function Qu{y, t). By using the recursion relation 
Eq. H it can be checked that Qu{y,t) satisfies the fol- 
lowing integral equation, 

Q„(2/,i) = G(y,t|0,0) 

-u f dt' r dy'G{y, t I y', t')V{y')Q^{y' , t'). (9) 

Jo J-00 

Using the definition of the propagator G, this integral 
equation can then be converted to a partial differential 
equation, 



dQu{y,t) _ ,2a^id^Qn{y,t) 



dt 



= af 



dy^ 



uV{y)Qu{y,t) (10) 



with Qu{y,t^ 0) = 6{y) and V{y) = sgn{y). 
We first make the scale transforms: 



y 



1 



z = — ; a = ut- Quiy, t) = —F{z, a). 

Substituting in Eq. Uu we get the following equation for 
F: 



dF d'^F OF 

a—- = a-—^ + az— \- \a — aV(z)]F 

da oz^ oz 



F{z,a = 0) 



V2 



2tt 



(11) 



where V{z) — sgn{z). This equation has the following 
series solution: 

_. 00 



1 °° 
F(z,a) = -=Vc„a"e-''Z?_„/„(z)e-^'/4 z > 0, (12) 
V 2tt „ 



where Dp{z) are parabolic cylinder functions. The coeffi- 
cients bn and c„ are to be determined from the boundary 
conditions, namely the continuity of both F and dF/dz 
at z = 0. The initial conditions determine &o = cq = 1. 
Using the boundary conditions we get: 



^ 6„a"e"i?_„/„(0) = Y. Cna"e-"^_„/,(0) 

00 00 

^ 6„a"e'^i?_„/„+i(0) = - ^ c„a"e-'^i^_„/„+i(0). (13) 



n=0 



n=0 



By expanding in powers of a and equating coefficients of 
all powers of a we finally obtain the following recursions 
for the coefficients fe„ and c„ 



E 



(-l)'"c„Z?_„./,(0) 



'p (n-m)! i?„„/„(0) 



(n odd) 



" ^ (-l)™c™i^_™/„+i(0) 

Z^ 7Z Z:XT~R TTTT^"-™ (" ^^^"^) 



„I^0 ("^™)' ^-n/a-n(0) 

-l)"c„ 



(14) 



where we have used few identities satisfied by the 
parabolic cylinder functions [Q. Now from Eq. ^ it 
follows that the moments ^k — {m^) satisfy. 



F{z,a)dz^Y. 



k=Q 



k\ 



-/ifc- 



(15) 



Finally, substituting the series solution for F(z, a) (Eq. 
|l2|) in the above equation we obtain: 



M« = n\ 



("l)'"c„ 



TT -"^-^ in — m)\ 

m=0 ^ ' 



in 



^-m/c-l(0) (16) 



for the even moments, while the odd ones vanish. The 
coefficents c,„ are determined from Eq. n4[ This thus 
gives an iterative scheme to generate all moments of the 
required distribution. 



III. AN ALTERNATE DERIVATION OF THE 
MOMENTS 

There is an alternate scheme to calculate the moments 
of the distribution P{m). This scheme assumes statis- 
tical independence of the successive zero crossing inter- 
vals of the process y{t) and was first used by Dornic and 
Godreche in the context of diffusion equation p3]. We 
however stress that while this assumption is only approxi- 
mate for non-Markov processes such as diffusion equation 
p4| , it is however exact for Markov processes such as the 
one we study in this paper. An additional complication 
in our case arises due to the fact that the average distance 
between zero-crossings vanishes. This is a standard result 
which is true for any Gaussian Markov process |32). In 
our calculations we introduce this average distance be- 
tween two consequtive zeros {1} as a cut-off parameter 
and then take the limit (l) ^ in the end. 

Consider a particular realization of the process y{t) 
ending at time t. Let at time t the process y have a posi- 
tive sign. Let i„ denote the time instant at which the nth 
zero-crossing takes place. Then the mean magnetization 
is given by: 

TO == -((t - tn) - {tn - tn-i) + ...) = 1-2^, where 



e 



1 r~ t 



Similarly if y{t) < then we get x = 2^ — 1. We note that 
at any t, the sign of y can be positive and negative with 
equal probability. Hence if we can find the distribution 
of ^, that of TO can be computed easily. Now in the 
logarithmic time variable T = log(t), we can write ^ in 
the form: 



= e-^(l-e-'"+e~'"-'"-i-f ...) 
= e^^Xn where 



(17) 



X„ = (l-e-'"+e-'"-'"-i +...), (18) 

/„ ~ T„ — T„_i and A is the time from the last zero- 
crossing to time t. The variables X„ satisfy Kesten re- 
cursion relations, 



Xn — 1 — e "X„_i. 



(19) 



One then assumes that the successive zero crossing inter- 
vals are statistically independent. In the long time limit 
the distribution of X is determined by the following set 
of equations: 

X = 77(1 -2C) + (l-r;) (2^-1) 
^ = e-^X 
X = 1-e-^X, 

where rj is an independent random variable that can take 
values and 1 with equal probabilities. Since one can 



compute the distributions of I and A, it is then straight- 
forward though tedious to compute all the moments of 



the mean magnetization, ^k 



recursively p4l 



As noted above, the mean distance between zero cross- 
ings, (/) vanishes and we introduce this as a cut-off pa- 
rameter. We now show that the /Jn are actually indepen- 
dent of (1). 

We first note that the Laplace transforms of the dis- 
tributions of I and A, which we denote by /(s) and g(s) 
respectively, are given by |24|: 



9(5) 



1 - (Og(5) 

1 + {l)g{s) 

Ms) 
■sil + {l)g{.s)y 



(20) 



where g{s) — s(l — sA{s))/2, and A{s) is defined 
as follows. Consider the normalized process Y = 
yit)/\/{y'^{t)). In the logarithmic time, T — log(t), this 
has a stationary autocorrelator, C{T =| Ti — T2 |) = 
g-"l^l. Now consider the autocorrelation function A{T) 
of the "signed'^ process, A{T) = {sgn{Y{0))sgn{Y{T))). 
The quantity A(s) is just the Laplace transforms of A(r). 
Using the fact that y{T) is Gaussian with a correla- 
tor C{T), the function A{T) can be easily computed, 
A{T) = (2/7r)sin-i[C(r)]. In our case, CiT) = e""!^! 
which finally gives. 



1 1 l^rS + aln 

s STT '-2a 2-' 



(21) 



where B[a, b] is the standard Beta function. 

It is now convenient to define the moments r„ =< 
X" > /(I -I- {l)gn)- Then taking the nth power of the 
equation X — 1 — e^'^X and using the expression for f{s) 
given in Eq. gO|, it can be shown that r„ 's are recursively 
generated through the following sets of equations: 

^/2n+l\, ,. 

2n-l 



2g2nr2n = - I] , (-1 



fc=l 



2n 



gkTk 



(22) 



with tq = 1. Note that all r„'s are independent of (l). 
The moments of ^ are then given by: 



< f ' >= q{n) < X" > = 



'-'9n^n 



(23) 



and are also independent of the cut-off (l). Finally 
the non-vanishing even moments of to can be obtained 
through: 



M2n - m - 1)'") 

and clearly do not depend on (l). 



(24) 



The final expressions of the first few even moments are 
as foUows ( see Eq. (B.4) in the appendix of ref. p4|), 



M2 = Ai 
/X4 = 1 - 



(l-3ii+4i2)(l-3i3 
I-2A2 



(25) 



We have checked that the moments /i„'s calculated re- 
cursively by this method are identical to those obtained 
by the first method in scction-II. 



IV. MOMENTS FOR SOME SPECIAL VALUES 

OF a 

In this section, we use the formalisms developed in the 
previous two sections to derive some explicit results for 
the moments of the distribution P{m). While the it- 
erative schemes developed in the previous sections are 
exact, it seems that for general a it is quite hard to ob- 
tain an exact closed form expression of /x„ for aribitrary 
n. They have to be determined only recursively. How- 
ever the equations simplify for some special values of the 
parameter a, for which not just the moments but the full 
distribution P{m) can be obtained explicitly. 

In order to see that the peak of the distribution shifts 
from m — ±1 for small a to to = for large a, it is 
natural to examine the two extreme limits a = and 
a — 00 for which fortunately we can obtain exact form 
of the distribution. Consider first a = 0. In this case it 
is somewhat easier to consider the second method used 
in section (III). It can be easily seen then that all the 
5„'s vanish while the moments r„'s remain finite. Thus 
from Eq. |2^ all moments of ^ vanish. Hence from Eq. |23 
we get /in = 1 for all even n. The same result can also 
be derived via the first method of section (II) by taking 
carefully the a ^ limit in Eqns. |l^ and |l6|. Immediate 
inspection then determines, 



-P(™) = -^ [S{fn - 1) + 5{m + 1)] 



(26) 



Now consider the other extreme limit, a = 00. 

= \/m\. 



14 that 



for a = 

In this case, one finds from Eq, 

Then the series in Eq. Hq just reduces to the expansion 

of (1 — 1)". Hence /x„ = for all n. This indicates that 

for a = 00 



P{m) = 5{m). 



(27) 



Another case where exact form of P{m) can be ob- 
tained is for a = 1/2. In this case, using the known values 
of the parabolic cylinder functions, it is easy to compute 
the first few terms of the series {c„,n = 0, 1,2, • • •} — 
{1,1,2,5,14,42,132,429,---} from Eq. |l|. We then 
make an ansatz, Cn = (2n)!/[n!(n -I- 1)!] and verify from 



Eq. |lj that it is indeed the solution for arbitrary n. 
Substituting this in Eq. |lj, we get 



M2n(a = 1/2) = 



(2n)! 
(n!)222»' 



(28) 



and the odd moments are identically zero. A little in- 
spection then shows that these are the moments of the 
distribution function. 



P(to) 



1 



ttVT 



(29) 



with TO varying in [—1,1]. We thus reproduce the well 
known p6| magnetization distribution for the ordinary 
random walk {a = 1/2). 

Unfortunately we were unable to get a closed form 
expression of /i„ for other values of a. For exam- 
ple, for a — 1/4, we get by solving Eq. |lj the 
first few terms of the sequence, {cn,n — 0, 1,2, • • •} = 
{1, 3, 72, 3663, 292824, 32227002, • • •}. We found however 
that this is not listed in the catalogue of known integer 
sequences |g5| and we could not guess any formula for 
this sequence. 

Thus as expected, the distributions of mean magneti- 
zation show a qualitative change in shape as a changes. 
As we go from small a to large a, the peak of the dis- 
tribution shifts from the edges to the center. This can 
be understood physically since for small a the noise be- 
comes small as time increases and the probability of zero 
crossing becomes negligible. On the other hand, for large 
a, the noise increases with time and the magnetization 
keeps changing sign and thus the most probable value 
gets peaked at to = 0. 

While obtaining an exact form of P{rn) is difficult for 
general a, there is no problem in obtaining the exact 
values of the moments of P(m) by using the recursion 
relations and the known values of the parabolic cylinder 
functions. In Fig. ^ we plot the moments for a = 1/4, 
1/2 and 3/4. 



V. CONCLUSION 

In this paper, we have studied the distribution of res- 
idence times or equivalently that of mean magnetization 
of a family of Gaussian markov processes parametrized by 
an index a which takes values continuously from to cx). 
We have shown that the shape of the distribution P{m) 
undergoes a qualitative change as a is increased from 
to CX3. For small a, P{rn) has peaks at the edges to, — ±1 
and has a minimum at ttt, = whereas for large a, the 
peak of the distribution shifts to ttt. = with minima at 
the edges to = ±1. This change in the ergodicity prop- 
erties of a stochatic process as one changes a parameter 
was first noted in ref. p5| in the context of diffusion equa- 
tion. The advantage of the process studied here, apart 



from representing various physical situations, is that the 
Markov nature of the process makes it possible to derive 
many exact analytical results. 

In this paper we have developed two alternate for- 
malisms to compute the moments of the residence times 
or mean magnetization. While both methods yield exact 
results for the moments, they do so only recursively. A 
closed form expression for the moments and hence that 
of the full distribution is possible only for some special 
values of the parameter a that characterizes the process. 
But unfortunately this special set of solvable values of a 
turn out to be the same for both these methods. Thus 
so far as the problem studied in this paper is concerned, 
both these methods are on equal footing. However there 
are other problems where the former method that gener- 
alizes Kac's formalism seems to have an advantage over 
the second method. We briefly mention below one such 
application. 

The general problem of a random walker in a space 
with moving boundaries has been well studied and has 
lot of applications |g^ . It would be interesting to study 
the residence time distribution in such problems. For 
example, consider a random walker moving in one di- 
mension and ask what is the distribution of the fraction 
of times spent by the walker in the region bounded by 
-|-oo and a point O that moves deterministically as xq (t) 
where xo{t) is some arbitrary function of i. For the spe- 
cial case when xo{t) = c\Jt where c is a constant, this 
problem can be solved by using the techniques presented 
in section-II of this paper. The calculations will be sim- 
ilar except that the potential Viz) — sgn(z) as used in 
Eq. O should be replaced by V{z) = 6{z — c). The cor- 
responding equations can be solved as before except that 
now the boundary conditions are to be applied at z = c. 
We note however that the second method illustrated in 
section-Ill does not seem to be easily generalizable to 
solve this problem. 

We conclude with one last remark. The magnetiza- 
tion distribution P{m) is a useful quantity to study for 
a generic stochastic process and contain in it many use- 
ful informations regarding ergodicity etc. However as is 
obvious from the efforts of this paper, exact analytical 
calculation of P{m) seems quite nontrivial even for the 
simple Gaussian Markov processes studied here. Thus at 
present, the only hope to compute P{in) for non-Markov 
processes which are richer and more abundant in nature, 
seems to be via numerical or approximate methods. 
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FIG. 1. In this figure the first few non- vanishing moments 
of P{m) are plotted, for a = 1/4, 1/2 and 3/4. 



